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ABSTRACT 

We discovered a large population of previously unknown Galactic H II regions 
by using the Green Bank Telescope to detect their hydrogen radio recombination 
line emission. Since recombination lines are optically thin at 3 cm wavelength, 
we can detect HII regions across the entire Galactic disk. Our targets were se- 
lected based on spatially coincident 24 /im and 21 cm continuum emission. For 
the Galactic zone —16° < I < 67° and \b\ < 1° we detected 602 discrete re- 
combination line components from 448 lines of sight, 95% of the sample targets, 
which more than doubles the number of known H II regions in this part of the 
Milky Way. We found 25 new first quadrant nebulae with negative LSR veloc- 
ities, placing them beyond the Solar orbit. Because we can detect all nebulae 
inside the Solar orbit that are ionized by O-stars, the Discovery Survey targets, 
when combined with existing H II region catalogs, give a more accurate census of 
Galactic H II regions and their properties. The distribution of H II regions across 
the Galactic disk shows strong, narrow (~ 1 kpc wide) peaks at Galactic radii 
of 4.3 and 6.0 kpc. The longitude- velocity distribution of HII regions now gives 
unambiguous evidence for Galactic structure, including the kinematic signatures 
of the radial peaks in the spatial distribution, a concentration of nebulae at the 
end of the Galactic Bar, and nebulae located on the kinematic locus of the 3 kpc 
Arm. 

Subject headings: Galaxy: structure — HII regions — radio lines: ISM — surveys 



1. INTRODUCTION 

HII regions are the formation sites of massive OB stars. Because the main sequence 
lifetimes of OB stars are < lOMyr, HII regions are zero age objects compared to the 
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age of the Milky Way. They thus reveal the locations of current Galactic star formation. 
They are the archetypical tracers of Galactic spiral structure. Their chemical abundances, 
moreover, indicate the present state of the interstellar medium (ISM) and reveal the elemental 
enrichment caused by the nuclear processing of many stellar generations. They provide 
unique and important probes of billions of years of Galactic chemical evolution (GCE). 
Individual Galactic H II regions are astrophysically important objects that reveal details of 
the impact of the star formation process on the ISM. Knowing the physical properties of H II 
region/photo-dissociation region/molecular cloud complexes provides important constraints 
on the physics of star formation and the evolution of the ISM. 

Moder n Galactic HII region surveys began with studies of the Palomar optical sur- 
vey plates (ISharplessI Il953l . Il959l ). Radio recombination lines (RRLs) from Galactic HII 
regions were discov e red in 1965 during the commissioning of the NRAO 140 Foot telescope; 
Hoglund fc Mezgerl (119651 ) detected H109a emission at 6 cm wavelength from M17 and 
OriA. Because the Galactic ISM is optically thin at centimeter wavelengths, RRL surveys 
were able to discover large numbers of H II regions distributed throughout the entire Galac- 
tic disk. Pioneering RRL surveys were done by, e.g., Reifenstein et al. 1970, Wilson et 
al. 1970, Downes et al. 1980, Caswell & Haynes 1987, and Lockman 1989. These sur- 
veys gave important insights into Galactic structure and the spatial pattern of massive star 
formation. Particularly noteworthy was the discovery of a metallicity gradient across the 
Galactic disk, made apparent by RRL measurements of HII region e lectron temperatures 
flWink. Wilson, fc Biegin3ll983l : lshaver et al.lll98~3l : buireza et al.ll2006h . T his Galactic-scale 
metallicity gradient placed strong constraints on GCE. By the time of the iLockmanl (119891 ) 
survey, however, almost all of the reasonably strong radio continuum sources, as revealed 
by contour maps often drawn by hand, had been observed; large angular scale surveys for 
discrete H II regions using cm- wave RRLs as tracers ceased. Despite these efforts, the cen- 
sus of Galactic HII regions was clearly incomplete. The advent of modern high-resolution, 
Galactic-scale infrared and radio surveys, coupled with the unprecedented spectral sensitiv- 
ity of the NRAO Green Bank Telescope (GBT), has made the GBT HII Region Discovery 
Survey (HRDS) possible. 



TARGET SAMPLE 



We assembled our target list from the following multi-frequency, large solid angle Galac- 
tic surveys: the infrared Spitzer Ga lactic Legacy Infr ared Mid-Plane Su rvey Extraordinaire 



(GLIMPSE: iBeniamin et all 120031 ) and MIPSGAL flCarev et al. 



Galactic Plane Survey made in 21 cm HI and continuum (VGPS: 



2009), the NRAO VLA 



Stil et al.ll2006h . the VLA 
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MAGPIS at 20 cm c ontinuum (IHelfand et al.ll2006h . and the NRAO VLA Sky Survey (NVSS: 
Condon et al.lll998l ). Our targets were selected by finding objects that showed spatially coin- 
cident 24 /iin Spitzer MIPSGAL and ~ 20 cm continuum emission, either from the VGPS or 
the NVSS. Our method was not automated, but instead relied on visual inspection of radio 
and IR emission maps. We removed known HII regions, planetary nebulae (PNe), luminous 
blue variables, and supernova remnants (SNRs) from the target sample using the SIMBAD 
database. 



This criterion is a strong indication that a target is emitting thermally (IHaslam fc Osborne 



19871 ; iBroadbent. Osborne. &: Haslamlll989l ). and therefore is likely an HII region or a PNe. 
Warm dust absorbs the far ultraviolet radiation from the exciting star(s) and re-emits in the 
IR, whereas plasma ionized by the same star(s) gives rise to free-free thermal emission at 
cm-wavelengths. It is certainly true that the radio continuum can result from a mix t ure o f 
free-free (thermal) and synchrotron (non-thermal) emission. iFiirst. Reich, fc Sofud (119871 ) 
showed that one can discriminate thermally from non-thermally emitting objects by using 
the ratio of the infrared to radio fluxes. The IR/radio flux ratio for HII regions is typically 
~ 100 times larger than that for non-thermally emitting SNRs, so it is easy to to differenti- 
ate between the two by eye. Our visual inspection of the IR and radio images should have 
eliminated all the SNRs from our target sample. Sources showing coincident mid-IR and 
cm-wave continuum emission almost invariably are thermally emitting: 95% of our sample 
targets show hydrogen RRL emission with line to continuum ratios of ~ 10 _1 which together 
suggests that our targets are emitting thermally. 



GREEN BANK TELESCOPE OBSERVATIONS 



Figure [TJ shows some representative GBT HRDS RRL detection spectra together with 
Spitzer MIPSGAL 24 /im images of the nebulae with contours of the VGPS 20 cm con- 
tinuum emission superimposed. (The MIPS detector is saturated in the central region of 
G32. 928+0. 607.) The sensitivity of the GBT and the power of its Autocorrelation Spec- 
trometer toge ther made the HRDS possible. To achieve high sensitivity, we used techniques 
pioneered by iBalserl (120061 ). who realized that there are 8 Una RRL transitions, H86a 
to H 93 a, that can be measured simultaneously by the GBT with the ACS at 3 cm wave- 
length (X-band). (The H86a t ransition, how ever, is spectrally compromised by confusing, 
higher order RRL transitions.) IBalserl (120061 ) showed that all these transitions can be av- 
eraged (after they are re-sampled onto the same velocity scale) to significantly improve the 
RRL signal-to-noise ratio, thus giving an extremely sensitive X-band Hn a average nebular 
spectrum. 
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This observing technique, coupled with the sensitivity afforded by the GBT's aperture, 
gave us unprecedented spectral sensitivity per unit observing time advantage compared with 
all previous cm- wavelength RRL surveys of Galactic HII regions. The vast majority of our 
detections took only a single OffOn total power observation. After Gaussian smoothing over 
5 channels to a resolution of 1.86 km s -1 chan" 1 (to be compared with the ~25kms _1 RRL 
typical line width), the r.m.s. noise for a single OffOn observation was typically ~lmJy. 
All details of the HRDS observing and analysis techniques, together with a catalog of the 
source properties, are discussed by Anderson et al. (2010, in preparation). 

We also used the Digital Continuum Receiver to measure the free-free continuum flux 
for all our HRDS targets. These show a power law distribution for fluxes > 120 mjy Be- 
low this the fitted power law overestimates the observed distribution which implies that 
the HRDS is complete to a flux l i mit o f > 120 mjy. Using the stellar fluxes given by 
Sternberg. Hoffmann, fc Pauldrachl (120031 ). we estimate that optically thin HII regions ion- 
ized by single 09 stars within the Solar orbit have flux densities > 120 mjy at 9 GHz. The 
HRDS thus should be finding all such HII regions and, furthermore, is capable of detecting 
03 stars at heliocentric distances of d Q ~20kpc ( Rubin Ill968l ; lAndersonll2009i ). 



4. DISCUSSION 

The GBT HRDS has doubled the number of known HII regions in the Galactic zone 
—16° < £ < 67° and \b\ < 1°. Each discovery gives a position and LSR velocity for 
the nebula. We detected 602 discrete Hydrogen RRL components from 448 target lines of 
sight, including ~65 infrared bubbles found in the Spitzer GLIMPSE survey. We found 
129 multiple velocity component sources (29% of the target sample). Figure [2] shows the 
longitude-velocity distribution of the HRDS nebulae, projected onto the Galactic plane. 

Although the HRDS found 448 HII regions, the number of physically distinct nebulae 
represented by the 602 discrete Hydrogen RRL components is not well known. Thermal radio 
sources can often be resolved into several apparently physically distinct emission regions, each 
having a somewhat different position and RRL velocity. Furthermore, the most massive star- 
forming complexes, W43 for example, are extended and fragmented into many sub-clumps 
of localized star formation, which together can ionize a very large zone. Many of our multiple 
velocity component targets may be detecting such low density ionized gas in addition to RRL 
e mission from anot her, physically distinct nebula. Here we follow the convention established 



by iLockmanl (119891 ) and assume that each of our 602 RRLs is produced by a distinct object. 
One should keep in mind, however, that the concept of a "discrete HII region" may not 
apply to many of the complex nebulae seen in the inner Galaxy. 
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Here we assume that all the HRDS sources are HII regions. Anderson et al. (2010, 
in preparation) show that the level of contamination in the HRDS sample by heretofore 
unknown SNRs, PNe, luminous blue variables, etc., is very small. This conclusion is based 
on considerations of galactic structure, scale height, RRL line widths, and the RRL line-to- 
continuum ratios (i.e. nebulae electron temperatures). 



Here we use the IPaladini et al.l (120031) cata log as a proxy for the sample of previously 
known H II regions. The IPaladini et al.l ( 120031 ) catalog of Galactic H II regions is a com- 
pilation of 24 Galactic HII region single-dish, medium resolution (~few arcmin) surveys 
covering the entire Galactic plane. There is, however, no definitive, complete compilation 
of all previously detected Galactic HII regions. Our HR DS nebulae are new, previously 
unknown HII regions because they are not listed in the IPaladini et al.l ( 120031 ) compila- 
tio n, the iLockman. Pisano. fc Howard! (119 96") sur vey, the Ultra- Com pact nebulae studied 
bv lArava et al.l J2002h . IWatson et ZlJ2003h . and ISewilo et al.l (120041 ) , nor are they in the 
SIMBAD database. 

The LSR veloci ty of each H RDS HII region maps into a unique Galactocentric ra- 
dius, Rq- Using the iBrandl (119861 ) rotation curve, Figure [3] shows the azimuthally averaged 
distribution of Galactic HII regions as a function of Rq- The filled histogram shows the 
H RDS nebulae; the op en histogram is the total HII region sample (1276 nebulae: HRDS 
& IPaladini et al.ll2003l ). The distribution of HII regions across the Galactic disk shows 



strong, narrow ( Rq < 1 kpc wide) peaks at Galactic radii of 4.25 and 6.00 kpc. There 



is still an overall dear th of HII regions within 4 kpc radius (jBurtonl Il976l ; iLockmanl 11981 
Bronfman et al.l 120001 ) . In the Rq = 2 to 4 kpc zone, however, the HRDS has found 75 



new nebulae, whereas IPaladini et al.l (120031 ) list only 56. The majority of these HRDS RRL 
components appear to be associated with the Galactic Bar at Rq ~ 3 kpc. The Galactic 
Bar and 3 Kpc Arm, however, produce large streaming motions that occur throughout the 
inner Galaxy, making all Rq deter minations that assume circular rotation uncertain inside 
4 kpc radius (IBurton fc Lisztlll993l ). 



Together, Figs. [2] and [3] show that the level of PNe contamination in the HRDS sample 
must be minimal. Because PNe are an old stellar population their Galactic orbits are well- 
mixed. PNe show, therefore, no structure in their Galactocentric radial distribution and 
their Galactic longitude-velocity distribution is a scatter plot constrained only by velocities 
permitted by Galactic rotation. Any PNe contamination of the HRDS sample must therefore 
be very small, otherwise these interlopers would suppress the unambiguous signal of Galactic 
structure seen in the HRDS Galactocentric radial and Longitude- Velocity distributions. 

Figure H] shows that the longitude- velocity distribution of the new HRDS and previously 
known HII region sample together now give unambiguous evidence for an ordered pattern 
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of Galactic structure. Doubling the census makes the contrast of the £-V features striking. 
The empty zones in Fig. H] are just as important in this regard as are the features seen. The 
sample of 1276 nebulae clearly shows the kinematic signatures of the radial peaks in the 
spatial distribution (Fig. [3]), a c oncentration of nebulae at the end of the Gala ctic Bar, at £ 
~30° and Vlsr ~+100kms _1 ([Benjamin et al.ll2005l ; IChurchwell et al.ll2009l ). and nebulae 
located on the kinematic locus of the 3 kpc Arm. 

The HII region Galactocentric radial distribution peaks at 4.3 kpc and 6.0 kpc have 
traditionally been associated with the Scutum-Centaurus and Sagittarius spiral arms. In 
the first Galactic quadr ant they impl y tangent point longitudes of 30° and 45°, respectively, 
for Rq =8.5 kpc. The iBrandl (119861 ) rotation curve terminal velocities for these directions 
are 107 km s -1 and 65kms _1 , respectively. The straight lines in Fig. [Hare the solid body loci 
defined by these tangent points and terminal velocities. These loci trace the over-densities 
in the HII region £— V distribution for £ > 0° quite well. These Northern tangent point 
longitudes show evidence for non-circular streaming motions of ~ lOkms^ 1 . 



Over 50 years after its discovery (Ivan Woerden. Rougoor. fc Oortl 119571 ) , the precise 
astrophysical nature of the M ilky Way's 3 Kpc Arm r emains enigmatic. The ellipse drawn in 
Fig. H] shows the locus of the ICohen fc Davied (119761 ) 3 Kpc Arm kinematic expanding ring 
model (Rg =4 kpc, V exp = 53kms _1 , V ro t = 210 km s -1 ). Discovered by the very first 21 cm 
HI surveys, the 3 Kpc Arm also contains considerable amounts of molecular gas (jBania 
19771 . Il980l . ll986h . Although the near side of the Arm, the segment between the Sun and the 



Galactic Center, is quite p rominent, the far side of the Arm was only recently discovered by 
Dame fc Thaddeud (120081 ). Because of this spatial symmetry and its extreme non-circular 
velocities, the 3 Kpc A rm provides strong evidence that the Milky Way is a m=2 barred 
spiral galaxy (jFux!ll999l ). 



Knowing that on-going star formation is occurring in the 3 Kpc Arm may help constrain 
theories of i ts dynamical origin. The 3 Kpc Arm lac ks large numbers of H II regions (ILockman 
19801 . Il98ll ). but there are some in it (jBanialll980l ). Although the HRDS found a few more 
nebulae (~ 10 more in the near side Arm, see Fig. [2]), it has not, however, found a substantial 
new population of HII regions along the 3 Kpc Arm l-V locus. The far side of the Arm, 
in the (£ < 0°, Vlsr > Okms -1 ) quadrant, remains almost devoid of nebulae along the Arm 
locus. Because the velocities predicted by the expanding ring model will blend with normal 
Galactic rotation velocities in the {£ > 0°, Vlsr > O kms" 1 ) and (£ < °, Vl sr < Okms -1 ) 
quadrants, the H II region clusterings seen along the I Cohen fc Daviesl (119761 ) ring locus are 
difficult to interpret. 

We were able to determine the distances to G38. 737-0. 140 and G48.551 -0.001 in Fig. [1] 
by using HI emission/absorption experiments to resolve the kinematic distance ambiguity 
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(As lAnderson fc Baniall2009l . did for the sample of previously known first Galactic quadrant 
HII regions.) Each nebula is at the far kinematic distance. The majority of our HRDS HII 
regions are unresolved with our 82" survey resolution. We are finding that many of these 
small angular diameter nebulae lie at the far kinematic distance. Earlier RRL surveys missed 
these nebulae because their weak continuum made them poor target choices. 

The HRDS found 25 first quadrant nebulae with negative LSR velocities. In the first 
Galactic longitude quadrant, a negative RRL LSR velocity unambiguously places the HII 
region beyond the Solar orbit, at large distances from the Sun, d & > 12kpc, in the outer 
Galactic disk, Rq > 9kpc. Prior to the GBT HRDS, there were only two negative-velocity 
H II regions known in the 18 ° < £ < 55 °, | b\ < 1 ° zone. The newly discovered Fig. 1 sources 
G38.651 + 0.087 at -40.0 km s^ 1 and G32.928 + 0.607 at -38.3 km s" 1 already match the size 
of the previous census. The t-V dist ribution of these negative veloci t y sou rces shows good 
agreement with 12 CO maps made by lDame. Hartmann. fc Thaddeud (120011 ). This region in 
£,K-space has traditionally been termed the "Outer Arm". 

Because of their location in the critical region beyond the Solar orbit at Rq ~ 9— 12kpc, 
these nebulae will provide new GCE constraints. The HRDS and follow-up GBT observations 
will allow us to derive the nebular electron temperature, T e , and helium abundances (Y = 
4 He/H). Because metals are the main coolants in the photo-ionized gas, both T e and Y are 
directly related to the distribution of heavy elements in the Milky Way. There are relatively 
few HII regions with accurately derived T e values, especially at the critical Rq ~ 10 kpc 
region. In the first Galactic quadrant, our 25 new HDRS nebulae can increase the T e sample 
size by a factor of 10. 

We detected RRL emission from 65 HII regions that are surrounded by mid-infrared 
bubbles. The Spitzer GLIMPSE survey found almost 600 obj ects in the inner Galaxy that 
have a ring-shaped morphology (IChurchwell et al.ll2006l . 120071 ). These are presumably bub- 
bles that are viewed in projection. These objects have simple morphologie s and "swept-up" 
neutr al material; they may be sites of triggered star formation (see, e.g., iDeharveng et al. 
2009h . IChurchwell et all J2006h argue that ~ 75% of GLIMPSE bubbles are caused by B- 
stars without detectable H II regions. Because of the large number of bubbles that we find 
enclosing HRDS sources, we speculate that nearly all GLIMPSE bubbles are caused by H II 



regions. 
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5. SUMMARY 

The advent of modern high-resolution, Galactic-scale infrared and radio surveys, e.g. 
Spitzer GLIMPSE/MIPSGAL and the VLA VGPS, coupled with the unprecedented spectral 
sensitivity of the NRAO Green Bank Telescope allowed us to make a major new discovery 
survey of Galactic HII regions. The GBT HRDS has doubled the number of known HII 
regions in the Galactic zone —16 < £ < 67° and \ b\ < 1 °. The census of HII regions, when 
enhanced by the HRDS, now shows in this Galactic zone a longitude-velocity distribution 
that gives unambiguous evidence for Galactic structure, including the kinematic signatures 
of peaks in the radial spatial distribution of nebulae, a concentration of nebulae at the end 
of the Galactic Bar, and nebulae located on the kinematic locus of the 3 kpc Arm. Doubling 
the Galactic H II region census makes the contrast of the t-V features striking. The empty 
zones in Fig. H]are just as important in this regard as are the concentrations of nebulae in 
l-V space. 

We found 25 new nebulae located beyond the Solar orbit, at large distances from the 
Sun, d & > 12 kpc, in the outer Galactic disk, Rq > 10 kpc. Because of their location, 
these nebulae will be important for future studies of the radial metallicity gradient in the 
Galaxy. Many of our new nebulae are seen as bubbles in Spitzer GLIMPSE images. We 
found 65 such objects and speculate that nearly all the Spitzer GLIMPSE mid-IR identified 
bubbles are H II regions. 

The HRDS nebular distances will be determined by usi ng HI emission /absorpti on ex- 



periments to resolve the kinematic distance ambiguity (see lAnderson &: Banial 120091 ). Be- 
cause we can detect all nebulae inside the Solar orbit that are ionized by O-stars, the GBT 
HRDS sources, when combined with existing H II region catalogs, will provide a more com- 
plete census of Galactic HII regions with known distances and physical properties, which 
is the fundamental database needed for ISM evolution studies of the Molecular Cloud/HH 
Region/Star Cluster/Supernova Bubble life-cycle. This will enhance our ability to study 
Galactic structure using spatial distributions and to constrain Galactic chemical evolution 
using spatial patterns of nebular metallicity. 



We thank those visionaries who came before us for the support to continue this sort 
of fundamental survey science. The National Radio Astronomy Observatory is a facility 
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Universities, Inc. LDA was partially supported by the NSF through GSSP awards 08-0030 
and 09-005 from the NRAO. 
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Fig. 1. — Left: GBT HRDS RRL spectra for four new HII regions. Shown are the average 
spectra for seven RRL transitions, H87o; to H93a, smoothed to 1.86 km s -1 resolution. 
Kinematic distances to these nebulae, d, were resolved using the RRL velocity and HI maps 
to study HI absorption of the HII region's continuum emission. The brightest HII regions 
show He and C recombination lines. Right: Spitzer 24 /im MIPSGAL infrared images for 
the same sources, together with contours of VGPS 20 cm continuum emission (1' resolution). 
Contours are drawn at 80%, 60%, and 40% of the peak emission. All images are 5' squares; 
each scale bar is 5 pc long. The GBT 82" (HPBW) beam is shown as a hatched circle. 
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Fig. 2. — Longitude- Velocity distribution of the 602 radio recombination lines found by the 
GBT HII Region Discovery Survey projected onto the Galactic plane. 
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Fig. 3. — Radial distribution of Galactic HII regions. The Galactocentric radiu s of each 
nebula, Rq, is derived from the observed RRL LSR velocity using the iBrandl (I1986I ) rotation 
curve. The filled histogram shows the 602 new nebulae found by the GBT HRDS. The 
open histogram show s the distribution of the 1,276 HII regions in the combined HRDS and 
Paladini et al.l (120031 ) samples. These histograms are averaged in Galactic azimuth, yet they 
show two significant, narrow ( Rq < 1 kpc) peaks at Rq = 4.25 kpc and 6.0 kpc. 
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Fig. 4. — Longitude- Velocity distribution of Galactic HII regions (the combined GBT HRDS 
& Paladini et al. 2003 samples) with kinematic Galactic structure loci. For clarity we have 
not plotted nebulae in the Nuclear Disk at large LSR velocities, \V\ > 150 km s -1 . Shown 
are the (£, v) loci for: (1) the circular rotation solid bodies for the peaks in the HII region 
Galactocentric radial distribution (straight lines: Fig. [3]); (2) the ICohen &: Davied (119761 ) 
3 Kpc A rm kin e matic ring model; and (3) the circular rotation terminal velocity derived 
from the iBrandl (119861 ) rotation curve. 



